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Abstract. In recent years tiie experimental study of neutrino oscillations has much 
contributed to our knowledge of particle physics by establishing non vanishing neutrino masses 
and by measuring or constraining the corresponding mixing angles. Within the domain of 
neutrino oscillations, the main goal of the next generation of facilities is the measurement of 
the mixing angle ^13 and the observation of leptonic CP violation, for which we have no hints 
at the moment. We discuss more in detail how various models of neutrino-nucleus cross section 
affect the forecasted precision measurement of 613 and the CP violating phase 5 [T]. 
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1. Summary of the charged current neutrino- nucleus cross sections 

Li the quasi-elastic regime (QE), the doubly-differential cross section, in which a neutrino 
carrying initial four- momentum k = {Ei^,\s.) scatters off a nuclear target to a state of four- 
momentum k = (-E'f,k ) can be written in Born approximation as follows: 



dndEi 167r2 |k| 

where Gp is the Fermi constant and Vud is the CKM matrix element coupling u and d quarks. The 
leptonic tensor is completely determined by lepton kinematics, whereas the nuclear tensor W^'^ , 
containing all the information on strong interactions dynamics, describes the response of the 
target nucleus, W^^'" = Y.X (0|^a^|^) {^\Ja\^) ^'^'^Hpo+Q-Px), where |0) and \X) are the initial 
and final hadronic states carrying four momenta po and px and is the nuclear electroweak 
current operator; the sum includes all hadronic final states. In the impulse approximation (lA) 
scheme, the nuclear current can be written as a sum of one-body currents, i.e. — t- J2i ^ 
while the final state reduces to the direct product of the hadronic state produced at the weak 
vertex (with momentum p ) and that describing the {A — l)-nucleon residual system, with 
momentum p?^: \X) — )■ |z,p ) (gi |7^, p-;^). 



1.1. The Spectral Function approach 

Following Ref. [2], the final expression of the hadronic tensor can be cast in the following form: 



= \ I d'pdEPip^E)--^ W^'^ip^q), (2) 



where = y|pp + m^ and the function P{p,E) is the target Spectral Function, i.e. the 
probability distribution of finding a nucleon with momentum p and removal energy E in the 
target nucleus. The quantity W^'^ is the tensor describing the weak interactions of the i-th 
nucleon in free space; the effect of nuclear binding of the struck nucleon is accounted for by 
the replacement = (i^, q) — t- g = (t'jq) with P = -Eip+qi — E^p^. The second argument in the 
hadronic tensor is p = {E^p^ , p). 

1.2. The Relativistic Fermi Gas 

The RFG [3j model, widely used in MonteCarlo simulations, provides the simplest form of 

the Spectral Function, PupcMip, E) = G{pf — p)S{Ep — Eb + E), where pp is the 

Fermi momentum and Eb is the average binding energy, introduced to account for nuclear 
binding. Thus, in this model pp and Eb are two parameters that are adjusted to reproduce the 
experimental data. For oxygen, the analysis of electron scattering data yields pf = 225 MeV 
and Eb = 25 MeV. 



1.3. The Relativistic Mean Field approach 

Within the Relativistic Mean Field approximation we refer to the model described in where, 
like in the previous cases, the nuclear current is written as a sum of single-nucleon currents. The 
wave functions for the target and the residual nuclei are described in terms of an independent- 
particle model. Then, the transition matrix elements can be cast in the following form: 

Jn(^, ?) = j dpTprip + q)JN{^, q)i^B{p) , (3) 

where ipB and ipp are the wave functions for initial bound and final outgoing nucleons, 
respectively, and is the relativistic current operator. In particular, the relativistic bound- 
state wave functions (for both initial and outgoing nucleons) are obtained as a solution of the 
Dirac equation, in the presence of the same relativistic nuclear mean field potential, derived 
from a Lagrangian containing a, lo and p mesons. 



1.4- The Random Phase Approximation 

The last model we want to take into account has been introduced in where the hadronic 
tensor is expressed in terms of the nuclear response functions treated in the Random Phase 
Approximation (RPA). The response functions are related to the imaginary part of the 
corresponding full polarization propagators and the introduction of the RPA approximation 
means that the polarization propagators are the solutions of integral equations involving the 
bare propagators and the effective interaction between particle-hole excitations. Within this 
formalism, the authors of [5j were able to show that multinucleon terms sizably increase the 
genuine charged current QE cross section in such a way to reproduce the MiniBooNE results [6] . 
The mechanism responsible for the enhancement that brings the theoretical cross section into 
agreement with the data is multi nucleon knock out, leading to two particle- two hole (2p2h) 
nuclear final states. In the following, we will refer to this "generalized" QE cross section as 
RPA-2p2h whereas we adopt the short RPA for the genuine QE cross section. 



1.5. Comparison of the cross sections 

To summarize this section, we present in Fig(T] a comparison of the five total QE cross sections 
for the ^^O fi~X process (left panel) and ^^O — )■ iJ,^X (right panel), in the energy range 
El, ~ [0,0.75] GeV. The curves have been computed using the dipole structure of the form 
factors and, in particular, a value of the axial mass close to mA_ ~ 1 GeV. As it can be easily 
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Figure 1. Total charged current QE cross sections for the i^jfO fi^X (left panel) and 
I'jfO — 7> (right panel) processes in the energy range Ey ~ [0,0.75] GeV. 

seen in the left panel, the RFC prediction sizably overestimates the SF, RMF and RPA results 
by roughly 15%, a fact that is well known to happen also for many other models with a more 
accurate description of the nuclear dynamics than the RFG approach. On the other hand, the 
inclusion of 2p2h contributions largely enhances the QE cross section in the RPA approximation 
for energies above ~ 0.5 GeV, although it is smaller than the RFG at smaller energies. For 
antineutrinos (right panel) the observed pattern is almost the same. 

2. The impact on the {9i^-5) measurement 

To estimate the impact of different models of the cross section on the measurement of ^13 and 
leptonic CP violation we choose a 7 = 100 /3-Beam facility as a representative example because 
the neutrino flux from such a facility spans up to ~ 0.7 GeV with the peak around 0.3 GeV and 
is thus mostly sensitive to the quasielastic region explored here. Details of the detector response 
can be found in Ref. [7j. It is important to notice that we are only using the quasielastic 
contribution to the neutrino cross section depicted in Fig(TJ As an illustration we have focused 
on the dependence on the nuclear model adopted of two different observables, namely the CP 
and 9i3 discovery potentials, defined as the values of the CP-violating phase 6c p and ^13 for 
which respectively the hypothesis of CP conservation 6c p = 0, ibvr or ^13 = can be excluded at 
3cr after marginalizing over all other parameters. The CP discovery potential is shown in the left 
panel of FigJJl where we superimposed the results obtained using the RFG cross section and the 
SF, RMF and RPA calculations. We clearly see that, for 6cp ^ ±90° (where the sensitivity is 
maximal) the RFG model gives a prediction which is around a factor 2 better than the SF, RMF 
and RPA models (which, as expected, behave almost in the same way) in sin^ 2^13 and around 
a 40% better than the RPA-2p2h. This is not surprising because the /3-Beam facility used in 
our simulations mainly probes energies smaller than 0.5 GeV, where the RFG is still larger than 
any other model (see Fig{T]). For the other points in the parameter space, the difference is less 
evident but still significant. The difference in the sensitivity to ^13 can be seen in the right panel 
of Fig. [2]where we show the results in the (sin^ 2^13, 5c'p)-plane. In this case the predictions of 
the SF, RMF and RPA models differ by up to a factor of ~ 60% compared to the RFG model 
for 6c p ^ ±90° while the difference is less pronounced for Scp ~ 0°. The RPA-2p2h results are 
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Figure 2. Left panel: CP discovery potential in the (^is, 5cp) plane. Solid lines refer to the 
SF model, dotted lines to the RFG, short- dashed lines to the RMF, dot- dashed to RPA and 
long-dashed to RPA-2p2h. Right panel: ^13 discovery potential. 

more similar to the RFG for 5cp ~ 0° and to the other models for 5cp ~ 180°. This is also 
evident in the right panel where we show the CP fraction. 

Finally, it is interesting to comment on the effect of using different nuclear models also in the 
simultaneous determination of ^13 and 5c p- We have performed such an analysis for the input 
value (^13, 5cp) = (0.9°, 30°) (not including degenerate solutions coming from our ignorance 
of the octant of ^23 [8J and in the neutrino mass ordering) finding that, using the RFG and 
RPA-2p2h models, we are able to reconstruct the true values of ^13 and 5c p within reasonable 
uncertainties, whereas with the other models we can only measure two distinct disconnected 
regions (the fake one around the value of ^13 and 5cp ~ 180°), which worsen the global sensitivity 
on those parameters. In conclusion, we have analyzed the impact of five different theoretical 
models for the neutrino-nucleus charged current QE cross sections on the forecasted sensitivity of 
the future neutrino facilities to the parameters ^13 and 5c p • We have found that the sensitivities 
computed with the FG are better, by up to a factor of 2, with respect to the other models. 
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